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Abstract 

To ascertain the recently observed boson at about 125 GeV by CMS and ATLAS Collaborations 
to be the Standard Model Higgs and determine its properties, there still needs more measure- 
ment on currently observed channels and more. In this work we reanalyze the processes of the 
Standard Model Higgs boson radiative decays to lepton pairs at one-loop level on the condi- 
tion of the LHC experiment. Result shows that the loop contributions for H — > = e, u) 
dominate over the Born level ones. In case the photon and leptons are hard, possessing ener- 
gies larger than 1 GeV, the branching fractions of H — > 11^(1 = e or fj.) processes are about 
two-thirds of the H — > u + ji~ process. Since the lepton-pair yields of the radiative processes 
mainly come from the Z-boson conversion, which may greatly suppress the background, we 
believe their signals should be observed in presently accumulated data or in the next run of 
the LHC experiment, provided the Standard Model Higgs is indeed light. 
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I. INTRODUCTION 



The recent discovery of the Standard Model(SM) Higgs-like boson around 125GeV 
by ATLAS [1] and CMS j^J Collaborations in Large Hadron Collider (LHC) experiment 
stirs the world with a great interest in high energy physics. Following, to ascertain the 
Higgs boson and further understand clearly its nature is one of the most important goals 
in high energy physics community To make sure the new finding is just the SM Higgs 
boson, but not others, more measurements on its decays are necessary. Therefore, to 
hunt for more LHC experiment accessible processes is also an urgent task for theorists. 

In ATLAS and CMS experiments, the SM Higgs-like particle was observed via its 
five decay processes, i.e., to 77, t + t and bb. All these decay channels 

are in two-body decay mode, or at least intermediately. In this paper, we propose a 
complementary channel, the three-body lepton-pair radiative decay processes H — > ii'y 
(£ stands for leptons), to reconstruct the SM Higgs-like neutral boson around 125GeV. 
Because of the internal heavy-quark and weak interaction gauge-boson contributions, 
these processes are not suppressed very much relative to the pure lepton-pair modes, 
and are found to be possibly observable within the datasets collected in 2011 and 2012. 

Before the LHC run, part of these radiative lepton-pair processes was analyzed [35], 
where the calculations were performed approximately and limited experimental environ- 
ment conditions were taken into account. In this work, we calculate completely the SM 
Higgs radiative decays to lepton pairs under the helicity basis, which may issue more 
information about the Higgs boson. Various physical cuts for reconstructing the Higgs 
boson in LHC experiment are applied in the calculation. The partial decay width with 
respect to the invariant mass of the muon pair is given in order to disentangle the SM 
Higgs boson with Higgs bosons in new physics models [6j. The H — > r + r~7 process is 
calculated beyond the Born level. Although the leading order tree process contributes 
dominantly, we find the one-loop process is very helpful in experimental measurement 
due to the Z-pole effect, which may greatly suppress the background. 
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Recently, Gastmans, Wu and Wu pointed out that the dimensional renormalization 
scheme may possess some shortcomings in the calculation of Higgs to two-photon decay 
process More explicitly, they believe that the W-loop contribution should decouple 

from the process while mathematically taking the zero-mass limit for W boson. However, 
hey found it is not the case in the calculation by dimensional renormalization scheme 
9], which has stimulated more discussions in the literature 10H14|. Since the triangle 
loop diagram in our calculation has a similar property to what in the Higgs to two 
photon process, we also check the W-loop decoupling theorem in this work. 

The remainder of the paper is arranged as follows. In Section II, the calculation 
procedure and the analytical results for the concerned processes are presented, and the 
W-boson decouple theorem is also confronted. In Section III, the numerical calcula- 
tion and the phenomenological study are performed. The last section is devoted to 
conclusions. 



II. THE DECAY WIDTHS CALCULATION 



A. The Tree Level Amplitudes 




FIG. 1: The Born level diagrams of Higgs boson radiative decay to lepton pairs. 



At Born level, there are two Feynman Diagrams for the Higgs boson radiative decays 
to lepton-pair processes, which are shown in Figure 1. For the convenience of following 
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discussion, the momenta of particles are assign as: p\ = pn, P2 = P-y, k\ = pe, and 
^2 = Pi- The Mandelstam invariants are s — (k± + k 2 ) 2 , t = (k 2 +P2) 2 , u = {k\ +P2) 2 , 
with s + t + u = 2mj + m 2 H . After some simplification, the tree level amplitudes are 
expressed as: 



_ e 2 mi (u(ki) $2 fo(k 2 ) - 2e ■ kiu(ki)v(k 2 ) 

■Mtree — ~ \ ^ 

2mwSw I 



+ 



u — m\ 

tt(fci) j> 2 fo(k 2 ) + 2e ■ k 2 u(k 1 )v(k 2 ) 
t — mj 



(1) 



B. The One-Loop Amplitudes 




(i) 





(2) 




FIG. 2: The typical loop diagrams of Higgs boson radiative decay to lepton pairs. 



For the concerned processes, the next-todeading order(NLO) contribution is very 
important for at least two reasons: first, the large couplings of heavy fermion and gauge 
boson to Higgs greatly enhance the decay rate in NLO. Secondly, the intermediate Z- 
boson decays may kick off a large amount of backgrounds and enable the experimental 



measurement more transparent. On the other hand, the Higgs radiative decay to lepton 
pairs itself is unique in determining the character of Higgs boson. 

The typical one-loop Feynman diagrams under consideration are shown in Figure 



2. In our calculation, both tree leve 
the Mathematica package FeynArts 



and NLO amplitudes are generated by virtue of 
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The analytical calculation of the amplitudes 



proceeds with the help of package FeynCalc [16]. Through out the calculation, the 
Feynman-'t Hooft gauge is adopted, and the one-loop Feynman diagrams are classified 
into four groups, i.e., 1) the group with triangle diagram and virtual photon directly 
coupling to ££; 2) group with triangle diagrams and virtual Z-boson directly coupling to 
££; 3) the group with W- and Z-box diagrams, and 4) the group with triangle diagrams 
but with photon emitted from one of the fermions. In practice, the contributions of loop 
diagrams with Higgs or Goldstone bosons directly coupling to leptons to the total decay 
widths are highly depressed relative to what from topologically the same diagrams with 
Z-boson coupling to the lepton pairs by a factor of which is much less than one and 
therefore can be neglected. 

In Feynman-'t Hooft gauge, there are 28 triangle diagrams in the first group. After 
some algebraic reduction and simplification the total amplitude of the processes in this 
group can be expressed as: 

M^M^ + M™ , (2) 

where the amplitudes for top quark and W boson(including Goldstone and ghost) in- 
duced loop diagrams read: 

H = TZ 5 ~ m Y N ( C 2 ^{(m 2 H -s)((m 2 H -s-Am 2 )x 

1 18n 2 m w sin 9 w s(m z H — s) 2 1 

C (m 2 H , 0, s, m 2 , m 2 , m 2 ) - 2) + 2s(B (s, m 2 , m 2 ) - B (m 2 H , m 2 , m 2 ))} x 
((m 2 H - s)u(h) fi(k 2 ) - 2e ■ Pl u(h) fcv(k 2 )) (3) 
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and 

AC = T7j 2 • f 7 2 ~ s)(m 2 w(™ 2 H - s)(6m 2 H - 12m 2 w - 5s)x 

lQTi z mw sm O w s{m H — s) z k 

Co(m^,0, s,m 2 ,m 2 ,m 2 ) - (m 2 H + Qm 2 v )(sB {m 2 H ,m 2 v ,m 2 v ) - sS (s,m^,m^) 
+ - s))u(A;i) ^(Afe) - 2{-2m 2 w {m 2 H - s)(-3m 2 w + 6m H + 4s)x 
C (m^, 0, a, m 2 , m 2 , m 2 ) - (m 2 , + 6m 2 w )(sB (m 2 H , m 2 w , m 2 w ) - sB (s, m 2 w , m 2 w ) 
+ m 2 H - s))e ■ piu(h) ^ 2 v(k 2 )} , (4) 

respectively. Since the bottom quark mass is about two orders smaller than the top 
quark mass, and the electric charge of bottom quark is half of top quark, it is reasonable 
to neglect the contribution of b-quark loop. The above amplitudes contain neither 
ultraviolet nor infrared divergences. The Cq and Bq functions are defined as in Ref. 

The value of the (^function and the finite part of the B function can be readily 
obtained by Looptools [121], and they have also been checked with analytical results. 

To check the W-loop decoupling theorem, we take a cut in the photon conversion 
diagram on virtual photon propagator, i.e., take the s to zero limit, and factor out the 
polarization vector, then the amplitudes turns to 

M %> = ~ 1fi^ m ,infl ro 2 « h + - rm%){2m 2 H + I2m 2 w + 12m 2 w (2m 2 w 

107T m\Y sm uy/ m H 

- m 2 H )C (m 2 H , 0,0, m 2 W} m 2 w } m 2 w )} . (5) 



When take the Co(m%, 0, 0, m 2 ^, m 2 ^, m^r) to be what given in Refs. [lOj, we get 
the same amplitude as obtained there. This again confirms that the W-boson loop 
contribution does not decouple from the total width in the zero mass limit. In fact, we 



ting 



ul, 



have checked that the W loop does not decouple in our interested processes too. Ta 
the s to zero limit for top-loop amplitudes, we get the same result as (4.11) in Ref. 
which in some sense provides one more check for our calculation. 

Similar to A^ 7 , there are also 28 triangle diagrams in group two with virtual Z boson 
directly coupling to the lepton pair. However, here the ultraviolet divergence appears, 
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which can be canceled out by adding the counter term in the on-shell renormalization 
scheme [3]. In the end, we get the renormalized and finite amplitude as: 

M z = M t z + M% , (6) 

where the amplitudes for top-quark and W-boson(including Goldstone and ghost) in- 
duced loop diagrams read: 

7Tl 2 C 4 N 

M * = ^rWfl cos^ m (m* C Twl ^TTf^ ^ 8sin2 ^~ 3 ){ 
o ( D7T sin u\\r cos uw m w \ m H ~ s) \s — m z + ll z m Z ) 

(m 2 H - s)((s - m 2 H + Am 2 )C (m 2 H , 0, s,m 2 ,m 2 ,m 2 ) + 2) + 2sB (m 2 H ,m 2 ,m 2 ) 

-2sB (s,m 2 t ,m 2 t )}((m 2 H - s)(w(fci) ^v{k 2 ) + (4 sin 2 6 W - l)u(kx) fifa)) 

-2e- Pl {u{h) i> 2 ^v{k 2 ) + (4sin 2 ^ - l)u{h) fov{k 2 ))) 

im 2 e 4 N c , 
967r 2 sin 3 9 W cos 2 Ow^iw ( m H ~ s )( s ~ m z + ^ z m z) 
(^ upa Pi P P2a((rn 2 H - s)C (m 2 Hl 0, s, m 2 , m 2 , m 2 t ) + 2B (m 2 H , m 2 , m 2 t ) 

- 2B {s,m 2 ,m 2 )))}e^(u(k 1 ) ll/1 5 v(k 2 ) + (4sin 2 ^ - l)u{k x )"( v v{k 2 )) (7) 
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and 

p 4 

j^w = f r 

1287r 2 m l y sin 3 8 W cos 2 9w{rn 2 H ~ s )( s ~ m | + ^ z m z) 

2(m 2 H — s)((3s — 5m^ + lSm 2 ^) cos 2 6^ + — 2s + m 2 ^) sin 2 x 
m 2 w C Q (m 2 H , 0, s, m 2 ^, m 2 ^, m 2 ^) + (cos 2 9 w (m 2 H + lOm 2 ^) 

- (m# + 2771^) sin 2 6 w )(m 2 H - s + sB (m 2 H , m 2 w , rr? w ) 

- sB Q {s,m 2 w ,m 2 w ))){u(k^ ^ 5 v{k 2 ) + {4sm 2 6 w - l)u(h) fi(k 2 )) 

+ , 

QAir 2 m w sin 3 W cos 2 6 w (m? H — s) 2 {s — m| + iT z mz) 

- 2(m 2 H — s)((6s — 5m^ + lOm 2 ^) cos 2 $w + (m 2 ^ — 2s — 2m 2 y ) sin 2 6W) x 
mfyCo^fy, 0, s, m 2 ^, m 2 ^, m 2 ^) - (cos 2 9 w (m 2 H + lOm^) 

- (m 2 H + 2m 2 w ) sin 2 6 w )(m 2 H - s + sB (m 2 H , m 2 Wl m 2 w ) 

- sB (s, m 2 w ,m 2 w ))}£ ■ piiuikx) ^ 2 n f 5 v(k 2 ) + (4 sin 2 9 W - l)u(fci) i> 2 v(k 2 )) 

647T 2 sin 3 Q w cos 2 ^(m 2 ^ — s)(s — m| + iT z m z ) 
2(m 2 H - s)(2cos 2 W - sin 2 w )C o (m 2 H , 0, s, m 2 w , m 2 w , m^) 
+ ^(s, m 2 w , mw) - B (m 2 H , m^, m 2 w ) }e ■ piu(ki)rfv(k 2 ) , (8) 

respectively. 

The calculation procedure for box diagrams is similar to what for triangle diagrams, 
except it becomes more tedious and complicate. However, the results for box diagrams 
are too lengthy to be presented here. Note that all box diagrams do not contain ultra- 
violet and infra divergences. 



C. Helicity amplitudes 



Notice that the helicity information on final states may tell more on parent particles, 
we perform also the calculation on the helicity basis for final states. Employing the 



helicity method given in Refs. 



19 



by introducing a light-like momentum q and a 



8 



space-like vector q% with constraints qo - qi — and qi ■ q\ — — 1, the helicity amplitudes 
can be then constructed as 

Mss'X = iV Tr[(# 2 - m,)(l - 7 5 ) j/ (#i + m*) A x ] , 

-M-s-s'A = AT Tr[(^ 2 - m/)(l + 7 5 ) ^o(#i + , 

-M- SS 'A = A^ Tr[(^ 2 - m,)(l - 7 5 ) ^ ^i(#i + , 

M S _ S , A = N Tr[{fc -m e ) ^{1- 7 5 ) ^ (#i + m £ )A x ] . (9) 

Here, N — 1/ y/l6(q ■ ki)(q ■ fc 2 ); ^4a = *A A ree +„4. A 5 ° p , which can be easily obtained after 
removing the spinors out from the above mentioned amplitude; s = 1/2 and s' = 1/2 
denote for the spin projections of I and I respectively; A = ±1 represents the helicities 
of the final photon. Thus there will be eight different kinds helicity amplitudes. 

Each polarization state can be calculated independently and the total unpolarized 
matrix element squared can be simply obtained by summing them up. I.e., 

\M\ 2 = \M SS , X \ 2 + \M^ S , X \ 2 + \M^ X \ 2 + \M S - S , X \ 2 . (10) 

In practice, we have checked that helicity method agrees with the traditional technique 
in getting the total unpolarized matrix element squared. 



III. NUMERICAL RESULTS 



A. Input parameters and selection cuts 



In numerical evaluation, the relevant inputs are taken be 2lJ, |22j : 

m t = 172.0GeV,m w = 80.39GeV, a(m z ) = 1/128, T z = 2.48GeV, 
m z = 91.18GeV,m„ = 0.105GeV ,m e = 0.51MeV,m T = 1.776GeV 
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In recently, the ATLAS and CMS Collaborations in LHC experiment observe a Higgs-like 
boson with mass of 126.0 ± 0.4 ± 0.4GeV |l| and 125.3 ± 0.4 ± 0.5GeV [2;], respectively. 
In this work we treat this boson as the SM Higgs and give it a mass of 125.5GeV. 
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TABLE I: Selection cuts taken in numerical evaluation. 





(ra«")cut 


(^| 7 )cut 


("if 7 )cut 


(E e ) mt (GeV) 


(E E ) cut (GeV) 


(£ 7 ) cut (GeV) 


cut I 


25ml 

(J, 


25ml 


P 


1 


1 


1 


cut II 


50m? 




50m?. 
p 


6 


6 


6 


cut III 


75mf, 


75ml 


P 


10 


10 


10 



Formally, the partial decay width can be expressed as 

dT = jk^ MU3du (12) 

Here, s, t and u follow the definition above Eq.(l). The phase space spans before 
imposing any cuts read as 

tmax = < _ ( ft _ ~2 _ <^)2 t = < _ ( /£ _ m 2 + ™H_Z2.f (13) 

and 

Smax = TTIh ' ^min = 4?7l£ . (14) 

To confront the theoretical calculation to the experiment measurement, certain phys- 
ical cuts should be imposed according to the condition of experimental facility. In 
LHC experiment, given that s > (m^) cut , t > (m| 7 ) cut , u > (mf 7 ) cut , E t > (E e ) cut , 
Ei > (Eg) cnt) and E 1 > (£/ 7 ) cut , where the last three inequalities are taken to be in the 
center of mass system of Higgs boson, three physical cuts are performed in our numerical 
calculation, as presented in Table [fl These cuts facilitate the experimental tagging of 
photon, electrons and muons, especially the cut III. For Higgs radiative decay to r pair, 
the lepton energy lower limits are taken to be the r mass in CUT I while keep the same 
as in Table Ufor CUT II and CUT III, and the invariant mass limits are 25m 2 ., 50m 2 , and 
75m 2 ., respectively, i.e., replace the m M by m T in Table [J The kinematic and energy cuts 
imposed can not only help to suppress the background, but also to distinguish leptons 
of the decay H — > 11 from leptons of the decay H — > 
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B. Numerical results and discussion 



With the analytic expressions calculated, the input parameters given and the physical 
cuts selected in preceding sections, the decay widths and branching fractions of a 125.5 
GeV SM Higgs boson radiative decays to lepton pairs can be readily obtained. The decay 
widths and branching fractions for various cuts are presented in Table HT1 For illustration, 
the invariant mass distributions of the decay mode H — > fx + fi~ / -f are presented in Figs. 
3 and 4, for combined and separated helicity cases, respectively. As for Higgs to e + e~7 
decay mode, the contributions from tree diagrams and the tree-one-loop interference 
term are negligible in comparison with the one-loop contribution, as also mentioned in 
Ref. {4]. Hence, the invariant mass distribution of H — > e + e _ 7 process tends to be as 
the NLO mass distribution of H — > fi + fi~ / y decay mode, that is the subtraction of the 
total contribution by the tree level contribution as shown in Fig. 3. For H — > r + r~7 
process, the dominant contribution comes from the leading order diagrams, due to the 
relatively large Higgs-r coupling. However, the subleading contribution for H — > T + r~--f 
process from the internal Z conversion to r pair tends to be meaningful with the help of 
Z-pole veto in experimental measurement. 

Numerical evaluation indicates that the contributions of those Box and triangle dia- 
grams, in which the photon emits from one of the fermions, are tiny in comparison with 
what from the first two groups. Furthermore, within the loop diagrams those with a 
W-loop triangle are primary for the decay processes. Since the triangle loop diagram in 
our calculation having a similar property to what in the Higgs to two photon process, we 
also check the W decoupling theorem by taking the zero mass limit, and find that W-loop 
contribution will not decouple from the process in the dimensional regularization. 

Result shows that the sum of the branching ratios of H — > /i + /i _ 7 and H — > e + e~7 
processes with cut I is about two times bigger than that of H — > ££££(£ — fi or e) 
processes j^J, where no any cuts at m# = 125 GeV is applied, and also lager than the 
branching ratio of H — > decay channel. Since H — > ££££(£ = \i or e) processes 

have already been observed in LHC experiment, the Higgs radiative decays to lepton 
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TABLE II: The decay widths and branching fractions with various cuts. Here, the mass of 
Higgs boson is taken to be 125.5GeV. 





r e - 


- e ~ 7 (10- 7 GeV)/ J Br(10" 4 ) 




> - 7 (10- 7 GeV)/Br(10- 4 ) 


r T+r - 


- 7 (10- 6 GeV)/£r(10- 3 ) 


cut I 




4.29/1.04 




5.67/1.38 




13.4/3.26 


cut II 




3.89/0.94 




4.57/1.11 




7.45/1.81 


cut III 




3.61/0.88 




4.09/0.99 




5.21/1.26 



TABLE III: The helicity decay widths with CUT I. Here, the mass of Higgs boson is taken to 
be 125.5 GeV. 





r e+e - 7 (10~ 7 GeV) /£r(10- 4 ) 




- 7 (10- 7 GeV)/£r(10- 4 ) 


r T+T - 7 (10- 6 GeV)/£r(10- 3 ) 


flip 


4.29/1.04 




4.30/1.04 


0.3/0.07 


non-flip 


0.00/0.00 




1.37/0.34 


13.1/3.19 



pairs are very likely to be also measurable. 

From Fig. 4 we notice that the contributions from muon and electron helicity non- 
flip processes are negligible in Higgs radiative decays. This is not strange, since in fact 
the helicity flip events are mainly induced by the loop diagrams. For tree diagrams, 
however, the largest contribution comes from the soft photon radiation where muon pair 
are produced back-to-back and photon is collinear to one of them, which leads to non-flip 
helicities of the muon pair. The detailed helicity results for decay width and branching 
ratio are given in Table IHIt where the mass of Higgs is taken to be 125.5 GeV. 

Furthermore, to enhance the tagging efficiency of final state particles and to depress 
the background, one can perform the measurement on muon pair with flip helicities and 
invariant mass in scope of 75 GeV to 100 GeV, where half of the CUT I events with 
photon energy ranging from about 20 GeV to 40 GeV will still be produced because 
of the Z-pole effect as shown in Figure 4. On the other hand, for this process the 
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FIG. 3: The \i\x invariant mass distributions of Higgs decay to fifi'y. Here, the premise of cut 
I is taken and the Higgs mass is set to be 125.5GeV. 

background seems to be understandable. The dominant background comes from the 
process qq — > fi + fi~j with photon emitted from muons and in soft region, which leads 
to flip helicities of muons. Most of them, however, are produced with a relatively large 
invariant mass and are out of the window. The same process but with photon emitted 
from initial quarks is relatively suppressed by factors of e\ or e 2 d . Hence this kind 
of background can be greatly depressed by measuring the muon helicity within the 
in scope of 75 GeV to 100 GeV of the muon-pair invariant mass. Another kind of 
background produced from those processes, like gg — > /i + /i~7, yeild muon pairs with 
flip helicities. However, generally these processes are higher order ones and involve loop 
diagrams, hence relatively suppressed by orders of strong coupling a s . Be that as it may, 
the background discrimination from the signal of the interested process will still be a 
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FIG. 4: The fip, invariant mass distributions of Higgs decay to with different muon pair 
helicities, the flipped ( — h and H — ) and non-fhpped( — and ++) cases. Here, the premise of 
CUT II is taken and the Higgs mass is set to be 125.5 GeV. 

challenge work. 

In comparison, processes H — > 7Z — > / -f£ + £~(£ = fi or e), where the Z boson is real, 
give the same final states as our concern and a branching fraction of 1.04 x 10~ 4 , half 
of what obtained in this work. In these processes, the emitting photon possesses a fixed 
energy of about 30 GeV and muons with energies varying from 30 GeV to 60 GeV in 
the Higgs center-of-mass system(CMS). 
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IV. CONCLUSIONS 



In this paper, the processes H — > tt^{t = e,fi,r) are investigated up to next-to- 
leading order in perturbative expansion in electroweak interaction. The branching ra- 
tions and decay widths are calculated in general and under helicity basis on condition 
of certain physical constraints, the requirements of minimum muon-pair invariant mass, 
etc. The CUTs may facilitate the experimental discrimination of the background from 
the signal. Numerical result shows that the branching ratio of H — > = e, fi) is at 

the order of 10~ 4 , about two times larger than that of H — > 'jZ — > , ~f£ + £~(£ — \i or e), 
while H — > TT7 process having a branching ratio of 10~ 3 . Our calculation indicates 
that although the Z-pole effect is overwhelming, the photon conversion to lepton pairs 
is not negligible, especially in case of a minimum lepton-pair opening angle. Since the 
Higgs radiative decays to electron and muon pairs have a branching ratio of order 10 -4 , 
by virtue of the helicity measurement and Z-pole constraint, it is expected to have a 
low background and the signals of these processes should be observed in presently accu- 
mulated data or in the next run of the LHC experiment, provided the Standard Model 
Higgs is indeed about 125 GeV. 

Recently, experimental progress has been achieved in measuring H — > rf process 
by CMS Collaboration [23(. Due to the large coupling of tau to Higgs boson, unlike 
the electron or muon pair radiative production, the tau pair radiative production in 
Higgs decay is mainly induced by the leading order contribution. However, the helicity 
measurement and Z-pole constraint may also be helpful to measure the Higgs radiative 
decay to tau-pair via loop contribution. 

Finally, we have also checked the W decoupling theorem in this work by taking the 
zero mass limit of the W boson, and find that the W-loop contribution does not decouple 
from the interested processes also in the dimensional regularization. 

Acknowledgements: This work was supported in part by the National Natural 
Science Foundation of China(NSFC) under the grants 10935012, 10821063 and 11175249. 



15 



[1] G. Aad, et al, the ATLAS Collaboration, Phys.Lett. B 716, 1(2012). 

[2] S. Chatrchyan, et al, the CMS Collaboration, Phys.Lett. B716, 30(2012). 

[3] A. Aabbasabadi, D. B. Chao, D. A. Dicus, and W. W. Repko, Phys. Rev. D55, 5647(1997). 

[4] A. Aabbasabadi, Wayne. W. Repko, Phys. Rev. D62, 054025 (2000). 

[5] Ana Firan and Ryszard Stroynowski, Phys. Rev. D76, 057301 (2007). 

[6] Chong-Sheng Li, Shou-Hua Zhu, and Cong-Feng Qiao, Phys. Rev. D57, 6928 (1998). 

[7] R. Gastmans, S. L. Wu, T. T. Wu, [arXiv:hep-ph/l 108.5322]. 

[8] R. Gastmans, S. L. Wu, T. T. Wu, [arXiv:hep-ph/1108.5872]. 

[9] G. Passarino and M. J. G. Veltman, Nucl. Phys. B160, 151 (1979). 

[10] William J. Marciano, Cen Zhang, Scott Willenbrock, Phys. Rev. D85, 013002 (2012) 

[11] Hua-Sheng Shao, Yu-Jie Zhang, Kuang-Ta Chao, JHEP 01, 053(2012) 

[12] Da Huang, Yong Tang, Yue-Liang Wu, Commun. Theor. Phys. 57, 427(2012) 

[13] M. Shifman, A. Vainshtein, M.B. Voloshin, V. Zakharov. Phys. Rev. D85, 013015 (2012) 

[14] Athanasios Dedes, Kristaq Suxho, [arXiv:hep-ph/1210.0141]. 

[15] T. Hahn, Comput. Phys. Commun. 140, 418 (2001). 

[16] R. Mertig, M. Bohm and A. Denner, Comp. Phys. Comm. 64, 345 (1991). 

[17] T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118, 153 (1999). 

[18] A. Denner, Fortschr. Phys. 41, 307(1993). 

[19] R. Kleiss and W. J. Stirling, Nucl. Phys. B262, 245(1985); Zhan Xu, Da-Hua Zhang and 

Lee Chang, Nucl.Phys. B291, 392 (1987). 

[20] Cong-Feng Qiao, Phys. Rev. D67, 097503 (2003); C.-H. Chang, J.-X. Wang and X.-G. 

Wu, Phys.Rev. D77, 014022 (2008). 

[21] K. Nakamura, et al, Particle Data Group, J. Phys. G37, 075021 (2010). 

[22] S. Dittmaier, et al, arXiv: 1101 .05931 [hep-ph], 

[23] Roger Wolf, Search for the SM Higgs Boson in Di-r Final States at CMS, Hadron Collider 
Physics Symposium, Kyoto, 2012. 



16 



